First and second sound in a highly elongated Fermi gas at unitarity 
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We consider a Fermi gas at unitary trapped by a higlily elongated harmonic potential and solve 
the equations of two fluid hydrodynamics at finite temperature. The propagation of sound waves 
as well as the discretized solutions in the presence of axial trapping are considered. The relevant 
thermodynamic functions entering the hydrodynamic equations are discussed in the superfluid and 
normal regimes in terms of universal scaling functions. Both first sound and second sound solutions 
are calculated as a function of temperature and the role of the superfiuid density is explicitly 
pointed out. The density fluctuations in the second sound wave are found to be signiflcantly large 
as a consequence of the coupling with first sound. Emphasis is also given to the comparison with 
recent experimental data. 

PACS numbers: 



I. INTRODUCTION 

Propagation of sound is one of the most exciting fea- 
tures exhibited by interacting many-body systems. It 
provides crucial information on the the dynamic regime 
obeyed by the system as well as on key thermodynamic 
quantities. The propagation of sound is particularly in- 
teresting in superfluids where two-fluid hydrodynamic 
theory predicts the occurrence of two different sounds 
P, 0: first sound, which is basically an iso-entropic wave 
and whose velocity is controlled by the adiabatic com- 
pressibility, and second sound, which corresponds to an 
entropy wave where the normal and superfluid compo- 
nents oscillate with opposite phase. Second sound has 
attracted much attention in the literature of superfluids 
mainly because the velocity of this peculiar sound is de- 
termined by the superfluid density. Actually in liquid 
^He the most accurate determination of the temperature 
dependence of the superfluid density is obtained through 
the measurement of the second sound velocity 0] . These 
data allowed Landau to establish the correct form of the 
spectrum of the elementary excitations of ^He, including 
the roton minimum Q . 

In ultracold atomic gases, below the critical tempera- 
ture for superfluidity, the propagation of sound has been 
the object of extensive theoretical and experimental work 
in the recent years. A peculiar feature of ultra-cold gases 
is that they are confined, the confinement being often of 
harmonic shape, which causes the discretization of sound 
waves in the form of collective oscillations. At the same 
time, in the case of highly elongated configurations, it 
is possible to investigate directly also the propagation of 
sound waves, by generating a perturbation in the center 
of the trap and subsequently investigating the time and 
spacepropagation of the resulting signal along the long 
axis dHl. 

Most of theoretical investigations in trapped atomic 
gases have been so far carried out at zero temperature, 
where only the first sound oscillations exist, both in Bose- 
Einstein condensed gas and in interacting Fermi gases. 



The agreement between theory and experiment is pretty 
good (for a review see, for example, Q and @). In 
the case of BEC gases these studies have permitted to 
check the validity of superfluid hydrodynamic theory at 
zero temperature [loj . In the case of interacting Fermi 
gases, along the BEC-BCS crossover, they have permit- 
ted to investigate fine details caused by the interactions 
and quantum statistical effects in the excitation spectrum 
of the discretized oscillations [ll|, including the beyond 
mean field Lee- Huang- Yang effect [Hj. The propagation 
of sound near zero temperature, and in particular the 
value of the sound velocity, have also been the object of 
systematic studies, and the general agreement between 
theory [l^ [l3| and experiments 0, Q is satisfying, both 
for Bose and Fermi superfluids. 

Much less is known about the behavior of sound at 
flnite temperature. This is due to various reasons. On 
one hand the realization of the collisional regime in the 
thermal component of a Bose gas, needed to apply the 
equations of hydrodynamics, is not obvious due to the 
very dilute and weakly interacting nature of the system. 
On the other hand the accurate control of temperature 
requires sophisticated experimental techniques. From the 
theoretical point of view the implementation of dynamic 
theories at finite temperature is much more difficult than 
at zero temperature jl5| , especially in the realistic case of 
trapped configurations. The situation becomes particu- 
larly challenging in the study of second sound, due to the 
intrinsic difficulties in generating and monitoring entropy 
waves and in providing accurate theoretical predictions 
for their behavior. 

Attempts to excite the relative motion of the conden- 
sate and thermal components of a harmonically trapped 
Bo se g as were first carried out in [l6| and more recently 
in [l3| where the use of denser samples permitted to ex- 
plore collisional damping effects. In Bose gases the space 
separation between the condensed and the thermal com- 
ponents caused by the trapping makes, however, the na- 
ture of second sound very different from the case of uni- 
form fluids where the two components fully overlap in 
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space. Furthermore, even in uniform configurations sec- 
ond sound in a Bose gas and in a strongly interacting 
superfluid, like ^ife behaves quite differently. In dilute 
Bose gases second sound actually reduces to the oscil- 
lation of the condensate [l^ , in agreement with Tisza's 
original idea. In a strongly interacting superfluid it in- 
stead corresponds to the oscillation of the gas of elemen- 
tary excitations relative to the superfluid component, and 
its observation consequently gives important information 
on the excitation spectrum. 

Differently from dilute Bose gases, Fermi gases at uni- 
tarity (where atoms occupying different spin states inter- 
act with a value of the scattering length much larger than 
the interparticle distance) behave like strongly interact- 
ing fluids and, in this respect, are more similar to liquid 
'^He. These systems actually exhibit novel and important 
features and are characterized by a universal thermody- 
namic behavior whose experimental determination has 
been the object of recent systematic studies [l^ [l^] and 
whose theoretical investigation represents a stimulating 
challenge of high interdisciplinary interest Q . Collisions 
in the unitary Fermi gas are more effective than in di- 
lute Bose gases, so that the non superfluid component of 
the system can easily achieve the hydrodynamic regime 
in a wide range of temperatures [2l| . Furthermore, due 
to their small compressibility, the superfluid and normal 
components strongly overlap in space also in the presence 
of external harmonic confinement. For the same reason 
the coupling between first and second sound is small [2^ , 
which permits to preserve the basic identification of sec- 
ond sound as an entropy wave, still giving rise to non 
negligible density fluctuations. 

The achievement of the hydrodynamic regime, as well 
as an efficient experimental excitation and investigation 
of sound waves is favored by the use of highly elongated 
traps. The discretized low frequency oscillations of first 
sound nature in these elongated configurations have been 
recently measured with high precision in the case of the 
unitary Fermi gas as a function of temperature [23j . ex- 
ploiting fine details of its thermodynamic behavior. The 
excitation and observation of wave packets propagating 
along the long (axial) direction has been also recently car- 
ried out at finite temperature 0, H^l ■ From the theoreti- 
cal side the use of highly elongated configurations allows 
for an important simplification of the formalism, through 
the formulation of the so called ID hydrodynamic ap- 
proach [25| whose derivation and implementation in the 
case of the unitary Fermi gas represent the main goal 
of the present paper. In this formalism the system pre- 
serves the applicability of the local density approximation 
along the axial as well as the radial directions, but ex- 
hibits a typical ID like behavior in the sense that the fluc- 
tuations of the relevant thermodynamic functions (such 
as the temperature, the velocity field and the chemical 
potential) depend uniquely on the axial coordinate [26j . 
This ID like behavior is the consequence of the effects 
of viscosity and thermal conductivity and is favored by 
a tight radial confinement. It is well suited to describe 



relevant experimental scenarios as we will discuss in the 
following sections. 

Our paper is organized as follows: In Sec. II we dis- 
cuss the main features of the ID hydrodynamic equations 
and the general conditions required for their applicabil- 
ity. In Sec. Ill we focus on the case of the unitary Fermi 
gas for which we derive the relevant 3D and ID ther- 
modynamic functions needed to solve the hydrodynamic 
equations. In Sec. IV we discuss the basic features of 
the variational approach which is implemented in Sec. V 
and VI, where we provide results for the first and second 
sound solutions, respectively. In particular in Sec. V 
we calculate the discretized frequencies which have been 
recently measured in [2^, while in Sec. VI we discuss 
the behavior of the second sound solutions both in the 
cylindrical geometry, where they take the form of sound 
waves, and in the case of axially trapped configurations, 
where the eigen-frequencies are discretized. In the same 
section we also discuss the coupling between the first and 
second sound modes. Finally in Sec. VII we draw our 
conclusions. 



II. TWO-FLUID HYDRODYNAMICS 

In this paper we consider atomic gases at finite tem- 
perature confined by a harmonic potential of the form 

T/ ^ 2 2,1 22 n\ 

and we will assume highly elongated configurations with 
trapping frequency satisfying the condition ujz ^ ^l, 
while m is the atomic mass. We will also consider the lim- 
iting Wz = case, corresponding to configurations with 
cylindrical geometry. Our aim is to discuss the dynamic 
behavior of a trapped Fermi gas of two spin species in- 
teracting with infinite scattering length (unitary limit), 
both below and above the superfiuid transition. We will 
discuss the propagation of sound in the cylindrical geom- 
etry {lOz = 0) as well as the discretized solutions with 
frequency of order lo^ ■ 

In [25j it was shown that, under suitable conditions of 
radial trapping, it is possible to derive simplified ID hy- 
drodynamic equations starting from the usual two-fluid 
Landau hydrodynamic equations defined in 3D [23| . The 
basic point for such a derivation is the requirement that 
both the velocity field along the long z-th axis and the 
temperature oscillations during the propagation of sound 
do not depend on the radial coordinates. This require- 
ment is justified in the case of tight radial confinement 
and is a direct consequence of the effects of viscosity and 
of thermal conductivity. The condition can be formulated 
in the simple form 

■q » mn„iUJ (2) 

where rj is the shear viscosity, n„i is the ID normal den- 
sity, obtained by radial integration of the 3D normal den- 
sity, and uj is the frequency of the relevant sound mode. 
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An analogous condition holds for the thermal conductiv- 
ity. In terms of the radial trapping frequency u}± the 
ID hydrodynamic condition ([2]) can be rewritten in the 
form uj w^T, where t is a typical coUisional time, here- 
after assumed, for sake of simplicity, to characterize both 
the effects of viscosity and of thermal conductivity [2^ . 
The above condition should be satisfied together with the 
usual hydrodynamic condition ujt ^ 1 [28| and is satis- 
fied by the very low frequency modes of the trapped gas, 
of order of uj^ , provided the radial trapping is sufficiently 
tight. It is worth noticing that the above ID condition 
implies that also the oscillations of the chemical poten- 
tial are independent of the radial variables. This follows 
psf from the condition of mechanical equilibrium along 
the radial direction dj_P + nd±Vext = and from the 
thermodynamic identity 

dP = sdT + nd^i (3) 

where s is the entropy density, n is the particle density 
and P is the pressure of the gas. Violation of the radial 
mechanical equilibrium condition would actually result 
in frequencies of the order aj_L rather than uiz- 

By radial integration of the 3D hydrodynamic equa- 
tions, and following the procedure described in [25j, one 
obtains the following ID hydrodynamic equations: 

mdtni+d,{i,) = Q (4) 



5tsi+a,(si0 =0 (5) 



mdtvl = -dMz) + VeAz)) (6) 



the system is large enough to safely carry out the radial 
integral using the local density approximation. 

The ingredients needed to solve the two fluid equations 
require the knowledge of the equation of state /i(r, n) and 
of the superfluid density. Theoretically, the calculation 
of the thermodynamic functions of the unitary Fermi gas 
is a great challenge due to the absence of a small coupling 
parameter. There are numerous efforts to develop strong- 
coupling many-body theories for such a sy stem (see [9| 
and |29| with references therein and [30l436{ ) . In order to 
probe the hydrodynamic behavior of the two fluid hydro- 
dynamic equations both the knowledge of thermodynam- 
ics below and above Tc is essential. Actually even below 
Tc the fluid is normal in the peripheral region where it 
approaches the classical regime. Due to the uncertain- 
ties of the theoretical calculations in relevant tempera- 
ture ranges, we have chosen the strategy of using, for 
the equation of state, the data available from the recent 
experimental analysis of the MIT team at unitarity [2^ . 
Universality can then be used to build the thermody- 
namic functions for all values of T and n (see Sec. III). 
Actually the experimental MIT data do not cover the 
whole range of temperatures and the information on the 
equation of state can be implemented and completed at 
high temperature by making use of the virial expansion 
( [29l[37| with references therein) and, at very low temper- 
ature, by calculating explicitly the phonon contribution 
which is known to give, in superfluids, the exact behav- 
ior as T — >• [l^. As far as the superfluid density is 
concerned, its present knowledge is rather poor and we 
will make use of simple phenomenological ad-hoc expres- 
sions in order to provide explicit predictions. Actually 
the investigation of second sound, which is particularly 
sensitive to the behavior of the superfluid density, is ex- 
pected to provide useful information and constraints on 
this quantity of fundamental interest. 



dtjz = -d,Pi-n^d,Ve.t{z) (7) 

where the terms ni, si. Pi are the radial integrals of their 
3D counterparts, namely the particle density, the entropy 
density and the local pressure, the integration accounting 
for the inhomogeneity caused by the radial component 
of the trapping potential ([1]). In the above equations 
jz = m{nniv^ + risiv^) is the current density, Usi and 
n„i are the superfluid and normal number densities re- 
spectively with ni = Usi + rini while and are the 
corresponding velocity fields. The continuity equation in 
Eq. Q expresses mass conservation. Equation ([5]) shows 
that the entropy of the fluid is carried by the normal fluid. 
In Eq. (UJ /ii = ii{T,n{r± = 0,z)) is the chemical cal- 
culated on the symmetry axis of the trapped gas and is 
determined by the equation of state of uniform matter. 
Eq. dni) fixes the law for the superfluid velocity, while Eq. 
([7]) is the ID version of the Euler equation for the cur- 
rent. Since we are interested only in the linear solutions, 
in the above equations we have omitted terms quadratic 
in the velocity. Here and in the following we assume that 



III. THERMODYNAMICS OF THE UNITARY 
FERMI GAS: FROM 3D TO ID 

In order to derive the relevant ID thermodynamic 
quantities needed to solve the hydrodynamic equations 
(|3][7|) let us flrst discuss the thermodynamic behavior of 
uniform matter. 



A. 3D thermodynamic functions 

At unitarity the s-wave scattering length diverges and, 
in uniform matter, the remaining length scales are the 
thermal wavelength 



At = \j2T:n^/mkBT (8) 

and the inter-particle distance n~-^/^. For the same rea- 
son the remaining energy scales arc flxcd by the temper- 
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ature T and by the Fermi temperature 



1 



(9) 



or, in alternative, by the chemical potential \i. It follows 
that at unitarity all the thermodynamic functions can 
be expressed [38| in terms of a universal function 
depending on the dimensionless parameter x = ^/kBT. 
This function can be defined in terms of the pressure of 
the gas, according to the relationship 



Using the thermodynamic relation n ~ 
density of the gas can then be written as 



(10) 

{dP/d^i)T: the 
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From Eq. ([TT|) one derives the useful expression 



T 



47r 



(11) 



(12) 



for the the ratio between the temperature and the Fermi 
temperature ([9|). 

In addition to the functions fp{x) and fn{x) it is also 
useful to define the function 



dx'fp{x'). 



which, as we will show soon, enters some relevant ID 
thermodynamic quantities. 

In terms of /„ and fp we can calculate directly the 
thermodynamic functions of the uniform Fermi gas at 
unitarity. For example, using the thermodynamic rela- 
tion S = V{dP/dT)^, we find 



S 



Nkf 



nks 2 /„ 



(14) 



while the specific heats at constant volume and pressure 
become 
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According to thermodynamics the ratio between Cp and 
Cv coincides with the ratio between the isothermal (kt) 
and the adiabatic (ks) depressibilities 
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(17) 



with kt and ks given, respectively, by : 



Kt 
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(18) 



(19) 



In Eqs. (|14j|16p we have introduced the entropy (s) and 
the specific heat [c^, cp) densities. 

The scaling function fp{x) (and hence the various ther- 
modynamic functions) can be determined through micro- 
scopic many-body calculations or extracted directly from 
experiments. In Fig. [l]and Fig. [2] we show, respectively, 
the equation of state ^i/kpT as a function of T/Tp and 
the universal functions fn{x) and fp{x) as a function of x, 
determined according to the procedures discussed in the 
following sections. In Fig. [3] we show the relevant ther- 
modynamic functions S/Nks, Cy/Nks and Gp/Nks as 
a function of T/Tp. 

Different from the previous thermodynamic quantities, 
the superfiuid density instead requires the knowledge 
of another independent function. According to dimen- 
sional arguments, at unitarity can be written as 



Arp 



(20) 



Its behavior is known at low temperature, in the phonon 
regime (see Sec. Ill C) [s^, and near the critical point 



(13) where it is predicted to vanish as Ug oc (1 — T/Tc 



c 



|2/3 

l40| . Here Tc is the critical temperature for superfluidity 
which, at unitarity can be written in the form 

Tc = aTp (21) 

with a a dimensionless universal parameter. 

B. ID thermodynamic functions 

Starting from the above 3D thermodynamic quanti- 
ties one can calculate the relevant ID quantities enter- 
ing the hydrodynamic equations (|4][7]), whose solution 
is the main goal of the present paper. In the pres- 
ence of radial harmonic trapping the chemical potential 
varies along the radial direction according to the law 
/^(''i) = Ml ~ (l/2)TOw^r^, predicted by the the local 
density approximation, so that one can easily reduce the 
radial integrals to integrals in the variable x. The follow- 
ing results hold for the ID density, pressure, entropy and 
specific heats per particle: 



ni{xi,T) 



dr±2TTr±n 



27r ksT 
nui'^ A|, 



fpixi), (22) 



Pi(xi,r)= [ drj_2nr^P=^^^^^fgixi) (23) 
J jnuj \ Arp 
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with 



Si(xi,T) 



= / dr±2nr±s ~ 



27r knTJ 



X 



fsi{xi) 



dxfs{x) 



(32) 
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Cpi(a;i 



dT 



4 fp{xi) 4 /„(a;i) 



(25) 



c„,(.05^(^ (26) 



where xi = fii/KsT is the value of the chemical poten- 
tial, in units of kgT, calculated on the symmetry axis of 
the trap, and Si = si/ni is the entropy per particle. 

Since Si depends only on the variable xi one finds that 
the adiabatic derivative of the ID pressure with respect 
to the ID density takes the form 



From the knowledge of the ID thermodynamic func- 
tions one can also easily calculate the equilibrium prop- 
erties in the presence of axial harmonic trapping, using 
the local density approximation ^ii{z) = hq — Vextiz) for 
the chemical potential along the z-th direction, with 
fixed by the normalization condition J dzni(z) = N. For 
example the ID density profile is available from Eq. ([2^ 
by replacing xi with [/io — {l/2)mojfz'^]/kBT. It is then 
natural to express the value of Xq = fj-o/ksT in terms of 
the Fermi temperature Tp""^ = {•iN)^/^huJho/kB of the 
3D trapped Fermi gas, where uiho is the geometrical aver- 
age of the three oscillator frequencies and N is the total 
number of atoms. One finds 



dPi 
dill 



5 ni 



(27) 



differently from the uniform case where one has, at uni- 
tarity, {dP/dn)s = i5/3)P/n with s = S/N. From the 
comparison between Eq. ([TU| and Eq. ([^^ one also finds 
the relationship 



(6/V^)-i/3( / dxixo ~ xy/^ux))-'/^ 

J — OO 

(33) 

This temperature scale will be used to discuss the tem- 
perature dependence of the discretized frequencies of the 
elementary excitations in the presence of ZD harmonic 
confinement. 
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P{r± = 0) 



(28) 



C. Phonon regime in the low T limit 



between the ID density and the pressure calculated at 
= 0. This relationship holds in the local density ap- 
proximation for a general fiuid radially confined with har- 
monic trapping [4l| . It actually follows directly from the 
radial integration of the general thermodynamic equation 
n = {dP/dfi)T. 

In Fig. |4] we show the relevant ID thermodynamic 
functions calculated function of the ratio 
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where 



ID 
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(29) 



(30) 



is the natural definition for the Fermi Temperature in 
ID cylindrical configurations [2^1 ■ T]P coincides with 
the usual 3D definition ^ of Fermi temperature if ni is 
calculated for an ideal Fermi gas at zero temperature. If 
one instead calculates ni in the unitary Fermi gas at zero 
temperature one finds the relationship T}P = ^^/^Tp 
where ^ is the Bcrtsch parameter (see next section). 

As concerns the ID superfluid density, starting from 
Eq. we find the expression 



At very low temperatures, corresponding to T << Tc, 
phonons provide the leading contribution to the ther- 
modynamic behavior of uniform superfluids. In this 
regime one can easily calculate the relevant thermody- 
namic functions introduced in the previous section. 

Starting from the expression [l^ 



Vn^kBT)' 
90(?ic)3 



(34) 



for the free energy associated with the phonon excita- 
tions in a uniform 3D superfluid, one can easily evaluate 
the other thermodynamic functions. In the above equa- 
tion Eq is the ground state energy, c is the T = value 
of the sound velocity, while V is the volume occupied 
by the gas. For the unitary Fermi gas one can write 
Eq ^ ^N^keTp and mc^ = Here ^ is the so 

called universal Bertsch parameter ([4^,0), accounting 
for the interaction effects of the unitary Fermi gas. Start- 
ing from p4p one can calculate the low temperature ex- 
pansion of the chemical potential /i — {dF/dN)T,v , pres- 
sure P = -{dF/dV)T.N and entropy S = -{dF/dT)v,N- 
One finds: 



nsiixi,T) 



dr±2TTr±ns 



2TT kuT 



fsiixi) (31) 



(35) 
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and 



NkB SO^Tf^ 



Using Landau criterion for the calculation of the 
phonon contribution to the normal density one also finds 
the result 



in order to determine the coefficient of the T^l"^ law the 
knowledge of the functions /„ , fp and fq are needed for all 
(36) values of x. This is physically due to the fact that in the 
radial integration the whole temperature range (and not 
only the large x phonon region) enters the calculation. 

It is finally interesting to calculate the low temperature 
expansion of the ID chemical potential. Equation ([25)1 
relates the ID density to the pressure of the gas calcu- 
lated on the symmetry axis so that the equation of state 
/ii (rii , T) corresponds to the equation of state of uniform 
matter as a function of P and T. We can consequently 
employ Eqs. p5ll36p to derive the low T expansion 



(37) 



(38) 



Using definition ((36)) one can also calculate the large x 
expansion of fp and hence of /„ . We find 



and 



5 37r2 



96 x' 



37r^ ^ 480 X 



(39) 



(40) 



The first terms in the expansions determine the T = 
value of the thermodynamic functions, while the second 
ones account for the first contribution due to the thermal 
excitation of phonons. 

An interesting consequence of the expansions pQMOp 
concerns the explicit low T behavior of the ID ther- 
modynamic functions. The ID entropy 1^^, as well as 
the ID specific heats (|25ll26p and the ID normal density 
nni = ni — Usi (see (Eq. I3ip ) exhibit a different T depen- 
dence as compared to the corresponding bulk quantities 
(see Eqs. P7ll38p ). In particular they behave like T^/^ 
as T — > 0, as can be shown by exploiting the leading be- 
havior at large x. For example the ID entropy behaves 



Si(T) 2TT{kBT) 



5/2 



(41) 



with 7 = dx'[^fp{x') - x' fn{x')\, the integral being 

convergent since {5/2)fp{x) — xfn{x) decays like x~^^'^ 
for large x 43|. Analogously the ID pressure and the ID 
normal density at low temperature behave as 



x/3^4 



f — )^ 



80^3/10^10 



(44) 



It is worth pointing out that, differently from the case 
of Eqs. ()41ll42p . the first contribution due to thermal ef- 
fects to III is determined by the phonon contribution and 
exhibits the typical dependence. It is also interesting 
to notice the opposite sign exhibited by the thermal cor- 
rection with respect to the bulk Eq. (j35|) which implies 
that these ID-like configurations will exhibit a different 
thermo-mechanical effect as compared to uniform gases 



D. Virial expansion at high T 

The high T behavior of the thermodynamic functions 
is determined by the virial expansion which takes the 
form of an expansion in terms of the fugacity (e^ ^ 1) 
and applies to very large and negative values of x. The 
hmction fp{x) can be expanded as: 



fp{x) = 2{bie^ + 62e2- + ...) 



(45) 



where bj arc the so-called virial coefficients and the fac- 
tor '2' comes from spin degeneracy. The value hi = 1 
is fixed by the classical equation of state, while theoret- 



ical calculations have provided the values 62 



3V2 



m 



and 63 = —0.29 [46[ for the second and third coefficients 
respectively. These values are consistent with the mea- 
surement of the equation of state at high temperature 

By taking the derivative of the pressure with respect 
to X, we obtain, for the function /n(x), the expansion: 



7muj , 



27rr 



^2) 



/„(a;) = 2(6ie- + 262e2--f ...) 



E. The intermediate temperature regime 



(46) 



and 



, , 27r(fcBr)5/2 m f^°" , . X , X 

nniiT) - ^ 2 ivT^) ^ / dx:y,,ix) (43) 
'mujj_ 2-1:1% J-oo 

with the quantity 1^,1(0;) vanishing as ^^[^-(f )^^^ at large 
x, in the phonon regime. The above equations reveal that 



Through high-precision measurements of the local 
compressibility, density, and pressure, the MIT team 
measured the universal thermodynamic behavior of the 
unitary Fermi gas with high accuracy both below and 
above the critical temperature for superfluidity overcom- 
ing, in particular, the problem of the direct measurement 
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of the temperature of the gas [20|. These measurements 
provide an important benchmark for many-body calcu- 
lations at finite temperature applied to this strongly in- 
teracting system. They have permitted, in particular, to 
identify the superfluid phase transition at the tempera- 
ture Tc ~ oTp with a = 0.167(13), corresponding to the 
value Xc = fic/ksT = 2.48. Also the relevant Bertsch 
parameter ^, giving the ground state energy in units of 
the ideal Fermi gas value, was determined with high accu- 
racy ( ^ = 0.376(4)). These values are in good agreement 
with the best theoretical predictions based on accurate 
many-body calculations (jl^, HH and [13] )■ Concerning 
the critical temperature it is worth stressing that in [20| 
its value was identified by exploring the peaked struc- 
ture exhibited by the specific heat at constant volume 
near the transition (see Fig. [3]), by taking explicitly into 
account finite resolution effects. This yields a value of 
Tc slightly higher than the value where the measured 
specific heat exhibits the peak. Figure [3] shows that the 
peak exhibited by the specific heat at constant pressure 
is even more pronounced, in agreement with the general 
behavior of the specific heats near a second order phase 
transition |47| . The matching between the values ex- 
tracted using the MIT data and the predictions provided 
by phonon thermodynamics of Sec. Ill A is also reason- 
ably good, especially as concerns the ID thermodynamic 
quantities. At high temperatures these experiments also 
confirm with high accuracy the validity of the virial ex- 
pansion (see Fig. [T] and Fig. [3]) so that, for the goals 
of the present paper, the equation of state of the uni- 
tary Fermi gas can be considered known with reasonably 
good accuracy at all temperatures. In the following we 
will adopt the MIT equation of state, together with the 
low and high temperature behavior discussed in Sec. Ill 
A and B, to implement the calculation of the frequency 
of the collective oscillations and of the sound velocities 
within the hydrodynamic formalism. 



IV. ID VARIATIONAL EQUATIONS 

The frequencies w corresponding to the solutions of 
the two-fluid hydrodynamic equations (|3]E]) with time 
dependence proportional to e~*'^* can be derived using 
the variational procedure 



where |4S 



Suj'^/Sun = 5uj^ /5Us = 



(47) 



( / dz[—{ — )rs,{5niY + 2ni( — )-s,5ni5si + 



rii dni 



dni 



nii^)niiSsiy])/ J dzm[nsiul + Hniul] (48) 

and Us, Un are the displacement field of the superfluid 
and the normal fluid, related to the velocity flelds via 
Us = Vs and u„ = u„. This variational procedure is the 
ID version of the 3D approach previously developed in 




FIG. 1: Equation of state (i/KbT versus T/Tp. The 
blue dash-dotted line corresponds to the phonon contribu- 
tion to thermodynamics (Sec. Ill C). The red filled cir- 
cles correspond to the experiment data in higher T regime, 
while the black solid line to the virial expansion in classical 
limit(Sec. Ill D). The green arrow indicates the critical point 
Tc/Tf = 0.167(13). The inset on the upper right corner is an 
amplification in the lower T regime. 




1 

X=|j/kgT 



2X3 



FIG. 2: Universal scaling functions /„ and fp as a function of 
the dimensionless variable fi/ksT. See Fig. [T]for the notation. 



45 1| . Keeping si constant in the derivatives of Eq. 
corresponds to considering ID isentropic transfor- 
mations. The density and entropy variations 6ni and 6si 
with respect to equilibrium are given, in terms of the dis- 
placement fields Un and Us, by Sni = —d^lrisiUs + riniUn] 
and (5s 1 = —UndzSi + (si/ni)dz[nsi{us — Un)]. Equations 
(|47ll48p hold in both uniform and trapped configurations, 
the effect of the trapping potential entering through the 
position dependent thermodynamic functions at equilib- 
rium. 

While in uniform configurations one can directly solve 
the hydrodynamic equations ([3117]), the use of the vari- 
ational procedure is particularly convenient in trapped 
configurations where analytic solutions of the full hydro- 
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FIG. 3: Entropy and specific heats per particle in uniform 
matter, evaluated using Eqs. (|14ll6p . In the lower panel, 
the red square-guided line corresponds to Cp/NkB; the black 
full-circle-guided line to Cv/Nks- 
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FIG. 4: The ID entropy and specific heats evaluated using 
Eqs. (|24I26|I . In the lower panel, the red dashed line corre- 
sponds to Cpi/Ais; the black dash-dotted line to Cvi/ks- 



dynamic equations are not available. 

In order to provide a first quantitative solution of tlie 
variational approach., we will make simple assumptions 
for the displacement fields of the two fluids. For first 
sound, which basically corresponds to a density wave, we 
will assume that the two fluids move in phase with equal 
velocity fields, i.e. Us = it,, = u. For second sound we 
will instead make the assumption that the total current, 
proportional to w„ri„i -l- Ms't-si, vanishes during the os- 
cillation. The coupling between the two modes will be 
discussed in detail in Sec. VI. Differently from the case 
of dilute BEC gases, where the coupling strongly affects 
the sound velocities in almost the whole temperature do- 
main [f^ , in the case of Fermi gases the coupling does 
not introduce important changes in the value of the col- 
lective frequencies, but is important because it releases 
the assumption that the total current vanishes in sec- 



ond sound, allowing for significant density fluctuations 
and hence providing imp ortant perspectives for its ex- 
perimental detection |52l |. 



V. FIRST SOUND 

Employing the first sound ansatz u„ = Ug = u the 
expression (j48p for the frequency to be used in the vari- 
ational calculation takes the simplified form 



(49) 



where we have employed the thermodynamic relation 
dzPi = -~nidzVext{z) holding at equilibrium (see Eq. 
([7])). In the unitary Fermi gas, where the ID thermo- 
dynamic relation (j27p holds, the variational procedure 
Soj^/5u = yields the following equation for the dis- 
placement field: 



7 y du 
— mojtz- 



dz 



7 Pi d^u 
5 ni dz'^ 



(50) 



The above equations explicitly reveal that the ID pres- 
sure Pi is the relevant thermodynamic quantity for de- 
scribing first sound dynamics. Equation (jSOp implies that 
the center of mass oscillation (dipole mode) characterized 
by the velocity field u = const is independent of the equa- 
tion of state and that the axial breathing mode [u oc z) 
takes the temperature independent value to ~ ^YljhbJz. 

For axially uniform configurations (w^ = 0) Eq. ([50)1 
predicts the propagation of sound waves with dispersion 
a; = cig and 



5 n\ 



(51) 



This differs from the sound velocity in uniform Fermi 
gases at unitarity, given by m(? = {^/'i)P/n, the differ- 
ence being caused by the presence of the radial trapping 
which gives rise to a different condition of adiabaticity. 
In Fig. [5] we show the value of the first sound velocity 
ci as a function of T/Tp^ using the thermodynamic re- 
sults for Pi/ 111 discussed in the previous sections. Using 
the expansion (|42|) for the ID pressure one finds that 
at T = the first sound velocity approaches the value 
cl = C^/^w|.ic/5 = Cw|/5 where v^i^ = 2mTlP and 
v'^ = 2mTp are, respectively, the 3D and ID Fermi ve- 
locities. The quenching of the sound velocity by the fac- 
tor •\/3/5 with respect to the bulk value S,vjp /3 was first 
pointed out in |14| . in analogy with a similar behavior 
exhibited by Bose-Einstein condensed gases [l^l ■ 

Let us now discuss the behavior of the discretized 
collective oscillations in the presence of axial harmonic 
trapping. At T = the ID pressure exhibits the po- 



sition dependence Pi/ni = (2/7) (/xq 



^muj^z' 



'), while 



at high temperatures, in the classical limit, one has 
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FIG. 5: ID first sound velocity in units of vjf calculated 
using Eq. (fSTj) . 



Pi/ni = ksT. Ill both cases, it is immediate to find 
that the sohitions of the hydrodynamic equation ()50p are 
polynomials of the form: = a^z^ -\- ak-2Z^~'^ + ... with 
integer values of k. At zero temperature one finds the 
following dispersion relation: 



bJt 5 



while in the high temperature limit one finds 



7fc + 5 



(52) 



(53) 



As expected, Eqs. ([52|) and (f53| coincide for fc = 
(center of mass oscillation) and k = \ (lowest axial 
breathing mode), while they predict different values for 
the higher nodal solutions. The result for the fc = mode 
follows from the universality of the center of mass oscilla- 
tion in the presence of harmonic trapping. The fact that 
the frequency of the lowest axial breathing oscillations 
does not depend on temperature is instead a peculiarity 
of the unitary Fermi gas. It is consistent with the exact 
scaling solutions exhibited by the two fluid hydrodynamic 
equations at unitarity (53| . The discussion then reveals 
that only the k = 2, 3... modes are useful in order to 
explore the effects of the temperature dependence of the 
equation of state. 

In order to provide a simple quantitative prediction 
for the temperature dependence of the k — 2 and k ~ 
3 frequencies we develop a variational approach to the 
solution of the hydrodynamic equations with the ansatz 
u = a2z'^+ao and u = a^z^ + aiz. This ansatz reproduces 
exactly the frequencies in both the T = and high T 
limits. Carrying out the variation with respect to the 
parameters characterizing the displacement fields, after 
a straightforward algebra one finds the results 



129*2 - 25 
45*2 - 25 



(54) 



440<3 - 252 
5(25*3 - 21)' 



(55) 



for the k — 2 and fc ~ 3 frequencies, respectively, where 
t2 = M0M4/MI and *3 = M2MQ/MI and we have intro- 
duced the dimensionless moments 



Ml 



dx{xo 



fn{x) 



(56) 



where xq is related to the value of T/Tp"^ by Eq. ([55]) . 
The integrals can be calculated using the data for the 
thermodynamic function fn{x) discussed in Sec. Ill 
which include the pro per interpolation between the ex- 
perimental data from [20[, the low temperature phonon 
regime as well as the classical regime, relevant to describe 
the low density region on the tails. 

The resulting predictions for the temperature depen- 
dence of the frequencies are shown in Fig. [HI together 
with the asymptotic zero temperature and classical val- 
ues as well as with the recent experimental results of [1^ . 
The results are plotted as a function of T/T^''"^ (see Eq. 
(|33l) ). The non monotonic temperature dependence in 
the higher temperature region is caused by the presence 
of the 2nd virial correction into the equation of state. 
We have also verified that the variational predictions for 
the collective frequencies of the fc = 2 and fc = 3 modes 
are practically indistinguishable from the exact numeri- 
cal solution of Eq. ([5(1)) . 

Using the same ansatz for the velocity field and the 
equation of continuity we can also calculate the density 
fluctuations of each mode given by 6ni ~ —dz[niu\ with 

^ - Ik^Th. ^ - 6k^Tu- equuibrium 
density profile ni{z) is available from Sec. Ill B, while 
the moments Mg are given by Eq. (|56| . In Fig. [7]we show 
the predictions for the density fluctuations calculated at 
T = for the fc = 1 and fc = 3 oscillations, while in 
Fig. [5] we show the comparison between the theoretical 
predictions for the k — 2 mode and the corresponding 
gerimental value at T = O.IT^™^ and T = 0.45Tj,™P 

The comparison between theory and experiment is in 
general quite satisfying, confirming the validity of the 
ID hydrodynamic approach used to predict the temper- 
ature dependence of these low frequency modes as well 
as the correctness of the equation of state employed in 
the calculation of the integrals M^. In particular the rel- 
atively small damping shown by experiments in the case 
of the fc = 2 mode confirms that the main assumption 
cj ^ ui\t needed to derive the ID hydrodynamic equa- 
tions is reasonably well satisfied. Preliminary experimen- 
tal data on the frequency of the fc = 3 mode also reveal 
a good agreement with theory, although in this case the 
observed damping is higher [54| . 
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FIG. 6: Frequency for the k = 2 (upper panel, Eq. ([54])) and 
k = 3 (lower panel, Eq. (|55|) ) first sound collective frequen- 
cies. Experiment data are from [S^]. The green lines are the 
theoretical predictions based on Eqs. (|54)l55p using the equa- 
tion of state of the unitary (solid) and ideal(dashed) Fermi 
gas. The thin horizontal dashed lines mark the zero-T su- 
perfluid limit (|52p and the classical hydrodynamic limit (|53[l . 
respectively. The red dash-dot vertical lines in (a) and (b) 
indicate the critical temperature. 
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FIG. 7: Density oscillations for the k = 1 (the blue dot-dashed 
line) and k — 3 (the red solid line) first sound modes at zero 
temperature. 



VI. SECOND SOUND 

Second sound corresponds to an out-of-phase oscilla- 
tion of the normal and superfluid components of the fluid. 
As a first ansatz we assume that the total current be 
zero {ji^z = "mriniv^ -\- mrisiv^ = 0) which implies that 
the oscillation is a pure temperature (or entropy) oscilla- 
tion, without any fluctuation of the density. Differently 
from first sound the superfluid density Ug plays a crucial 
role in the propagation of second sound. In superfluid 
'^He the measurement of the second sound velocity has 
actually provided the most accurate determination of Ug 
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FIG. 8: Density oscillations for the k — 2 collective oscil- 
lation (upper figure) and equilibriurn profiles (lower figure) 
at two different temperatures, from [23|. In this figure the 
Fermi temperature corresponds to the definition T^""^ = 
(3Ny^'^hu}ho/kB introduced in the text. 



as a function of temperature Q. At present the knowl- 
edge of Us is rather poor in the unitary Fermi gas and we 
will make use of simple ansatz for a first estimate of the 
frequency of the second sound oscillations. Under the as- 
sumption that the total current vanishes, the expression 
(|48l) for the frequency of second sound to be used in the 
variational calculation takes the simplified form: 



UJ 



J rir, 1 ^ 



(57) 



and the variational condition SuP' /5us = yields the fol- 
lowing equation for the displacement field of the super- 
fluid component: 



UJ^Us 



Si 



^ dT d_ 



mn\ dz ds 



)] (58) 



The above equations reveal that the key thermodynamic 
quantities characterizing the propagation of second sound 
are the ID density, entropy, specific heat and super- 
fluid density. The presence of axial trapping is indirectly 
present through the value of the equilibrium quantities. 

From Eq. (j58p one immediately recovers the second 
sound velocity for an axially uniform system by consid- 
ering a plane wave solution of the form e"^^ for Us . One 
flnds u! = C2q with 



2 rpSi "-si 

TOC2 = J 



(59) 



For uniform configurations it is actually possible to 
solve exactly the two-fluid hydrodynamic equations (|4][7]) 
and it is interesting to check the accuracy of the approx- 
imate prediction (|59p . Using straightforward thermody- 
namic relations it is actually possible to show that the 
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solutions for the sound velocity emerging from the HD 
Eqs. (|3][7]) with Vpxt{z) = should satisfy the following 
equation 



2rl .3Pi 
C — U— 51 

m ani 



1 UsiTsl^ , 1 UsiTsl dPi 



m nniCyi 



(BO) 



yielding two solutions for the sound velocity, correspond- 
ing to the first (ci) and second (C2) sound velocities. An 
accurate expression for the lower solution (second sound 
velocity) is derived under the condition 



C2 Cp\ 



Cvl 



Cvl 



< 1 



(61) 



In this case one finds 



T 



si^ nsi 



(62) 



which differs from Eq. ((59)) because of the presence of 
the specific heat at constant pressure rather than at con- 
stant density. The two specific heats actually exhibit a 
different behavior for temperatures close to Tc (see Fig. 
3]). When T — >■ the specific heat at constant pressure 
and at constant volume coincide and, as a consequence 
of the temperature dependence of the ID thermodynamic 
functions in the low temperature regime (see Sec. Ill B), 
the second sound velocity vanishes like VT, differently 
from what happens in the bulk where it approaches the 
value ci/-\/3 |18l |. At finite temperature expression (|62p 
is very accurate in reproducing the lower solution of (|60p 
for all temperatures, the condition ((6T|) being always well 
satisfied. The above discussion then reveals that second 
sound can be regarded as an oscillating wave at constant 
pressure, rather than at constant density, as previously 
assumed in the derivation of ((59)) . This is the conse- 
quence of the coupling with the first sound mode. Differ- 
ently from the second sound velocity, the velocity of first 
sound is instead negligibly affected by the coupling, pro- 
vided condition ((6ip is satisfied and is consequently very 
accurately described by Eq. ((5T)) for all temperatures. 

The coupling between first and second sound has the 
important consequence that the density fluctuations dur- 
ing the propagation of second sound are not vanishing. 
Actually, under the condition C2 ^ ci , the ratio between 
the relative density and temperature fluctuations takes 
the simple expression 



Sni/ni 
ST/T 



n\ 01 



7 fnfq 

2lf 



(63) 



55|. The ratio §^ is 



and turns out to be negative 
shown in Fig. [9] and should be compared with the result 
^"^^1^ = ■^ir§^)si = I characterizing the propagation 
of iirst sound, where the derivative is calculated at con- 
stant entropy. It is remarkable that the ratio ((S5)) is sig- 
nificantly large in a useful range of temperatures, thereby 
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FIG. 9: Ratio (|63p between the relative density and temper- 
ature fluctuations calculated for the ID second sound in the 
axially uniform configuration. 



opening realistic perspectives for the observation of sec- 
ond sound waves, by looking at the density fluctuations 
during the propagation of sound. 

In Fig. (TOjwe show the prediction for the second sound 
velocity ((5^ using two different models for Ug. A first 
model employs the formula n^/n = (1 — T jTc)'^!'^ , cor- 
rectly accounting for the behavior near the critical point. 
A second model employs the phenomcnological expres- 
sion Tis/rt = 1 — (T/Tc)* which exhibits, at low temper- 
ature the correct T"* behavior, although the coefficient of 
the law is about 8 times larger than the one predicted 
by the phonon contribution to the normal density (see 
Eq. ((38)) '). As expected the second sound velocity de- 
pends in a crutial way on the choice of the model for Us 
so that the measurement of the C2 is expected to provide 
unique information on its temperature dependence. 

In the presence of harmonic trapping along the z-th di- 
rection, the second sound modes are instead discretized 
and we use the variational approach (|57p to obtain first 
estimates for the collective frequencies. The lowest fre- 
quency mode of second sound nature is expected to be of 
dipolar nature for which we make the simplifying assump- 
tion that the velocity field of the superfluid component 
is constant in space. In Fig. [11] we show the resulting 
prediction for the temperature dependence of the lowest 
second sound mode, using the two models for the super- 
fluid density discussed above. We have checked that the 
inclusion of higher order terms in the polynomial ansatz 
for Us introduces minor corrections (less than 10 %). 

An important feature of the second sound frequency 
is that it vanishes when the temperature approaches the 
critical value. This result differs from the one predicted 
in 3D isotropic configurations jsTi . i56 | and can be under- 
stood noticing that an estimate for the discretized fre- 
quency can be obtained using the expression lo ^ Ciq with 
q ~ where Rs^z is the the size of the superfluid 

along the z-th direction. On the other hand the main 
temperature dependence of the second sound velocity, as 
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FIG. 10: ID second sound velocity in the axially uniform con- 
figuration, calculated using Eq. (|62p and two different models 
for the superfiuid density: the red dashed-line corresponds to 
the phenomenological ansatz: ris/n = 1 — (T/Tc)* while the 
black dash-dotted-line to the choice: Uajn = (1 — T jTcf'l^ . 

T Tc, is given by the the ID superfluid velocity that 
behaves like n^i ~ UgR"^ ^ and is hence proportional to 
the bulk superfluid density calculated in the center of the 
trap and the square of the size of the superfluid along the 
radial direction. Since both and the ratio ±/ Rs,z 
vanish as one approaches the transition temperature, the 
second sound frequencies will vanish too. 

Differently from the case of the uniform ID configura- 
tions in the presence of axial trapping the coupling be- 
tween first and second sound discretized solutions cannot 
be easily calculated in an exact way. It can be never- 
theless estimated using a variational approach. To this 
purpose we will look for solutions of the variational hy- 
drodynamic equations in the form Ug = au^^^ + u^s'^ and 
Un = au*-^) for the superfluid (us) and normal 

velocity fields, respectively. Here u^^-* = Un^ = u'p cor- 
responds to the velocity field of the first sound solutions 

(2) (2) 

discussed in Sec. V, while Un and ul arc the velocity 
fields of the uncoupled second sound solutions discussed 
above. By inserting the ansatz in (pS)) we find, after a 
straightforward calculation, the expression 

„2, ,2 I , ,2 K2 „ Ul,2 

= ' (64) 

for the collective frequency as a function of the vari- 
ational parameter a, where K\ = (1/2) J dzmniUi, 

K2 = (1/2) / fizm(ui^^)2ninsi/n„i and Ui,2 = 
(2/5)T/ dz§|L 9£i!iii|i!V!i:ii ^nd we have used the iden- 
tity 'T-i(^-)si = |T holding at unitarity. By imposing 
the variational condition Slj'^ /5a = we find the result 

uj = (65) 



FIG. 11: Frequency for the lowest discretized second sound 
mode in an axially trapped configuration with uJz <^ oj±. 



for the frcc^ucncy of the two coupled modes. When ap- 
plied to uniform matter, using result ([51]) and (|59|) for 
the uncoupled first (wi) and second (102) sound frequen- 
cies, the above procedure reproduces exactly the two de- 
coupled solutions given by roots of Eq. As an 
example of application in the presence of harmonic trap- 
ping we have considered the coupling between the dipole 
second sound solution discussed above (corresponding to 
ui^-* = 1 and Un^ = —uf^nsi/uni and the k = 2 first 
sound solution discussed in Sec. V. The k = 2 mode is 
actually the lowest frequency first sound mode that can 
be coupled to the dipole second sound mode. The nu- 
merical calculation shows that the changes in the value 
of the second sound frequency caused by the coupling are 
very small. 

Let us finally point out that the discretized second 
sound oscillations discussed above are expected to be 
more damped than the first sound ones discussed in the 
previous section. The reason is that the thermal con- 
ductivity in the normal phase tends to infinity near the 
transition point [STf and is consequently large near the 
boundary between the superfluid and the normal phases. 
This is expected to result in the penetration of the tem- 
perature oscillation into the normal phase, resulting in 
an increase of damping. 



VII. CONCLUSIONS 

We have provided a systematic discussion of the two 
fluid hydrodynamic behavior exhibited by the unitary 
Fermi gas in the presence of a highly elongated harmonic 
potential. The main achievements contained in the paper 
are summarize below. 

i) We have presented an exhaustive discussion of the 
relevant 3D and ID thermodynamic functions, like the 
pressure, the entropy and the specific heats at constant 
density and at constant pressure, whose knowledge is rc- 
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quired in order to solve the hydrodynamic equations. The 
thermodynamic functions are identified using the most 
recent experimental data obtained at MIT, through the 
introduction of universal scaling functions which empha- 
size the universality of the unitary Fermi gas. The match- 
ing of the MIT data with the low T behavior of the 3D 
thermodynamic functions fixed by the thermal excita- 
tion of phonons and with the high T virial expansion has 
been explicitly discussed. Particularly interesting results 
concern the behavior of the ID quantities which are cal- 
culated by radial integration of the 3D thermodynamic 
functions using the Local Density Approximation. The 
behavior of the ID thermodynamic functions at low tem- 
perature is not uniquely fixed by the thermal excitation 
of phonons as happens in uniform superfluids, but in- 
volves also the thermal regimes at higher temperature. 
Their temperature dependence at low T has been explic- 
itly calculated. 

ii) We have solved the ID hydrodynamic equations de- 
rived in j25j using a variational formulation of the hydro- 
dynamic equations. Explicit results are given for both 
the first and second sound modes. While the first sound 
solutions are basically determined by the ID adiabatic 
compressibility the second sound solutions are sensitive, 
in addition to the entropy and the specific heat, to the 
superfluid density of the system, a rather elusive quan- 
tity which cannot be determined by the knowledge of the 
equation of state of the system and whose knowledge in 
strongly interacting Fermi gases is still rather poor. 

iii) We have provided results for both ID uniform and 
axially trapped configurations. In the first case the so- 
lutions of the hydrodynamic equations take the form of 
sound wave whose velocity has been systematically in- 
vestigated for both first and second sound, employing 
different models for the superfluid density. In the pres- 
ence of axial trapping the lowest excitations take the 



form of collective oscillations whose discretized frequen- 
cies are calculated function of temperature. The 
theoretical predictions for the flrst sound solutions are 
compared with recent experiments carried out both be- 
low and above the critical temperature for superfluidity. 

iv) An important feature emerging from our studies is 
that in highly elongated configurations the coupling be- 
tween first and second sound makes the latter mode an 
oscillation at constant ID pressure, rather than at con- 
stant ID density and an explicit formula for the result- 
ing density fluctuations has been derived as a function of 
temperature. This has the important consequence that, 
except at very low temperature, the density fluctuations 
characterizing second sound are sizable, thereby making 
this mode observable in experiments. 

Open questions to address in future works concern the 
damping of the collective modes caused by viscosity and 
thermal conductivity and a more quantitative check of 
the applicability of the ID hydrodynamic description em- 
ployed in the present paper. The conditions of applica- 
bility of this approach, i.e. the independence of the fluc- 
tuations of temperature and of the velocity field on the 
radial coordinate, are actually ensured by the crucial role 
played by viscosity, thermal conductivity in the presence 
of tight radial trapping. 
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